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Abstract. There is now strong evidence that the current energy density of the Universe is dominated by dark 
energy with an equation of state w < —1/3, which is causing accelerated expansion. The build-up of structure 
within such Universes is subject to significant ongoing study, particularly through the spherical collapse model. 
This paper aims to review and consolidate progress for cosmologies in which the dark energy component remains 
homogeneous on the scales of the structures being modelled. The equations presented are designed to allow for 
dark energy with a general time- varying equation of state w(a). In addition to reviewing previous work, a number 
of new results are introduced: A new fitting formula for the linear growth factor in constant w cosmologies is 
given. A new formalism for determining the critical density for collapse is presented based on the initial curvature 
of the perturbation. The commonly used approximation to the critical density for collapse based on the linear 
growth factor is discussed for a general dark energy equation of state. Virialisation within such cosmologies is 
also considered, and the standard assumption that energy is conserved between turn-around and virialisation is 
questioned and limiting possiblities are presented. 
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1. Introduction 

Analysis of the distance-redshift relation using high red- 
shift Type la supernovae has led to the discovery that 
. the expansion of the Universe is curren tly accelerating 

■ ijRiess et al.lll998l IPerlmutter et al.lll999j) . This suggests 
. that the dominant contribution to the present-day energy 

■ budget is a component with equation of state w < —1/3, 
called "dark-energy" . Combining measure ments of CMB 
fluctuations (dp Bcrn a rdis et al . 2000; Lee et a ll |200ll: 
Halverson et al.ll2002t IPearson et alJ l2003t IScott et all 



20031 iHinshaw et all kool " with measurements of the 
clustering of present day galax ies |Percival ct al. 200ll 



Percival et al. JP02; 



energy (e.g. 


3fstathiou et al.' 


'2002! 


Socrgcl ct al. 


2003; 


Tegmark ct al. 



0.3, with the remaining contribution made up of dark en- 
ergy. 

The nature of the dark energy is the source of much 
debate. Perhaps the most straightforward candidate is a 
positive cosmological constant A with equation of state 
parameter w = —1. This simple picture forms a special 
case in a broader class of models where the dark en- 
ergy is the manifestation of a scalar field slowly rolling 



down its potential. In the limit of a completely fl at po- 
tential, these models lead to w = — 1 (WcttcricQ Il988l: 
IPeebles fc RatrJlQSStFRatra fc Pee'bielll988^b A subclass 
of quintessence r nodels with constant —1 < w < —1/ 3 
was proposed bv lCaldwell. Dave fc Steinhardtl l)l998albl) . 
Many other forms have been proposed for the shape of 
this potential, leading to an equation o f state parameter 
that is dependent on the scale factor fsee lPeebles fc Ratral 
I2OO3I for a review) . 

Given the plethora of models, many parameterisations 
have been introduced i n order to observationally constrain 
w{a) (see IJohril l2004l for a review). However, there is 
obviously a limitation imposed by any parameterisation 
l|Bassett. Corasaniti fc K unz 20 0^ and, in the present pa- 
per, we have tried to be as general as possible in the equa- 
tions presented and discussed. In order to demonstrate the 
effects of the dark energy equation of state in numerical 
example s, we pay particular attention to the p arameteri- 
sation of Ijassal. Baela fc PadmanabhanI l)2004() 



w{a) ^ Wo + wia{l — a). 



(1) 



The present day value of wja^i = wq, and its derivative 
dw{a)/da\a=i = —wi will serve to demonstrate the pos- 
sible cosmological signatures of a wider class of models. 
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We also consider models that assume 'w{a) = w is con- 
stant. Qualitatively, many of the effects of general w{a) 
models can be predicted by interpolating between models 
with constant w. In fact, for many models, by the look- 
back time by which w{a) has evolved significantly from its 
present day value, the cosmological significan ce of w(a) is 
great ly reduced (although see, for example, iMaor et all 
120021) . Obviously, this also signifies that an evolving w{a) 
is harder to obser vationally constrain than constant w 
l)Kinat et 3,1.1120021) . If the dark energy equation of state 
only varies slowly with time, then observational predic- 
tions are well approximate d by treating w{a) — w as a. 
constant with 



J daftxia)w{a) 
J da Q,x (a) 



(2) 



critical overdensity for collapse at present day (Section^, 
and as a function of time (Section Isj. This critical over- 
density is then used to determine the mass function 
(Sectional) and the rate of structure growth (Section [Tj). 
Finally we consider the virialisation of perturbations, and 
the difficulties associated with these calculations in dark 
energy cosmologies (Section ISJ. 

2. dynamics of cosmological models 

It is assumed that the dark energy has an equation of 
state relating its pressure px and density px given by 
Px = w{<i)px- For general w{a), the dynamical expansion 
of the Universe is specified by the Friedmann equation 



^^^'^) O -3 I o -2 I o /(a) 
— -rpy— — ilMa +\IkCL +\lxa \ 



Hi 



(3) 



where ^x is the dark energy density relative to the c ritical 
density (see also lDave. Caldwell fc Steinhardtll2002l) . 

The equation of state for the dark energy does not 
uniquely define the behaviour of this component. The 
formation of structure is also dependent on the sound 
speed of the dark energy which limits its clustering 
properties. In the original formalisrn for quintessence 
ijCaldwell. Dave fc Steinhard3ll998albi) . the dark energy 
component has a high sound speed which means that it 
can cluster on the largest scales, but does not cluster on 
the scales of galaxy clusters and below. 

Consequently, the dark energy only affects the mat- ^or constant w, /(a) = -3(1 + w). For the 

ter power s pectrum llMa et aljUooi. and th e CMB parameterisation w{a] = + w^ia(l - «) of 

anisotropics (ICaldwell. Dave fc Steinhardill998bl) on very l -Jassal, Bagla fc Padmanabhan| 12004D, 



where 51^- = (1 — SIm — f^x) is the curvature con- 
stant, iJ(a) = a/a is the Hubble parameter with 
present day value i/o- /(a) is calculated by solving the 
conservation of energy eq uation for the dark energy 
d{pxa^ )lda = -ipxa^ ijCaldwell. Dave fc Steinhardd 
Il998bl) . giving px oc a^'^°'\ where 



/(«) = 

In a 



In a 



[l + w{a')]d\TLa' . 



(4) 



large scales. For the special case of a cosmological con- 
stant, w = — 1, the clustering of the dark energy is not 
an issue as the energy density in perturbations always re- 
mains at the background level. Obviously, for w ^ —\, 
the clustering properties of the dark energy strongly af- 
fect the build-up of structure in the Universe. In partic- 
ular for spherical perturbations, the linear growth rate, 
critical overdensity for collapse, and details of subsequent 
behaviour and virialisation are all dependent on this prop- 
erty. In this paper we follow the majority of current lit- 
erature and only consider a non-clustering dark energy 
component. However, we do note that models in which 
the dark energy clusters on small scales are beiiig dis- 
cussed with iiicreasing frequency iIHu fc ScrantonI 12004 
Nunes fc Mota '2004': H annestadI l2005t iManera fc Motal 



/(a) = -3(l + u;o) + 



21na 



{l-af 



(5) 



The evolution of the matter density ^Mio) and dark 
energy density fix (a) are given by 



Q.M{a) 



fix (a) 



(6) 



Nunes & Mota 


'2004: 


2005: 


Maor & Lahav. 



In general, for the variables used in this paper, if no 
dependence is quoted for a given quantity (e.g. flAi), it 
should be assumed to be calculated at present day. If in- 
stead explicit dependence is given (e.g. flMia)), the quan- 
tity is assumed to vary with epoch. The exception to this 
rule is w{a), where if no dependence on a is given, we 
additionally assume that w{a) — w is constant in time. 

We start by reviewing the behaviour of different cos- 
mological models, and the critical parameters that sep- 
arate models with different properties (Section [SJ as this 
can be related to the behaviour of spherical perturbations. 
We then calculate the linear growth factor (Section|3J), the 



It is immediately apparent that the equation of state 
of the dark energy has a strong affect on the behaviour of 
these equations, even for models where w{a) is constant. 
For example, for w = —1/3, the dark energy terms in 
Eq.|21cancel, leaving a cosmological model with an expan- 
sion time history that behaves exactly as an open Universe 
with matter density f^Af, alt hough the s pace-time geom- 
etry differs (e.g. review bv IPeebles fc R.atra 2003i) . For 
w < —1/3, at late times the dark energy will dominate 
the dynamics of the expansion, with the epoch of transi- 
tion from matter to dark energy domination dependent on 
w. Decreasing the equation of state w from —1/3 smoothly 
interpolates between the open Universe model, cosmolog- 
ical constant (A model) and extrapolates beyond. 

The asymptotic behaviour (either forwards or back- 
wards in time) of different cosmological models is delin- 
eated by models with so-called critical parameters. These 
separate cosmologies that predict future recoUapse, ex- 
pansion forever and models that do not start at a big- 
bang, but instead "loiter" at early times, asymptotically 
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tending towards a stable solution such as Einstein's static 
model for A cosmologies. For these critical models, the 
zero points of (da/dt)'^ = E'^{a)a'^ = should occur at 
turning points of this equation. For constant w, differen- 
tiating Eq. 13 we see that the critical cosmologies should 
have E'^(a)a^ = when ijChiba. Takahashi fc Sugivamal 
120051) 



M 



(7) 



Substituting this value of the scale factor back into Eq. O 
gives an equation for the critical parameters. 

For w = —1, this re duces to the well known cubic 
equation 



27 
T 



(8) 



that can be solved analytically ijFelten fc Isaacmarll986() . 
A detailed discu ssion of the behaviour of A models is given 
bv lMolesI lll99j^ . 

For general w{a), it is straightforward to numerically 
determine the critical parameters from Eq. |21 and its 
derivative. The different regions in (f^M^^x) space are 
shown in Fig. ^ for wq = —4/3,— 1,-2/3 and Wi — 

— 1/3, 0, 1/3 using the parameterisation of Eq. ^ Critical 
values of ^Im and fix at a = 1 are shown by the solid 
black lines. 

For comparison, in Fig. ^ we also plot various cosmo- 
logical tracks (grey lines). For general w{a) models, the 
cosmology is not uniquely specified by flMi^i) and fix (a), 
as we also need to specify a. The cosmological tracks in 
Fig. ^for fix > have a — 1 when flM = fix- For gen- 
eral models, the tracks can therefore cross, and can move 
into regions of (flM{a-),flx{a)) space that would have been 
excluded at a = 1. For the panels where wi ^ 0, the cos- 
mological tracks should simply be though of as examples 
of possible tracks that go through particular flM{a) and 
fix (a). 

For wi = 0, Wo < — 1 models favour a larger region of 
parameter space for which there is no big-bang, while wq > 

— 1 leads to recoUapse for models in a larger region in {flM, 
fix) space, reflecting the fact that a smaller amount of 
matter is required to collapse the Universe before the dark 
energy dominates. For wq > —I cosmologies, a future sin- 
gularity is predicted where the scale factor and Hubble pa- 
rameter dive rge, leading to a "big-rip" rather than expan - 
sion forever dCaldwell. Kamion kowski fc Weinberdl2003() . 
For furt her discussion of the behaviour of wi =0 cosmolo- 
gies, ^ ee Ide Arauiol ll2005l) : IChiba.' Takahashi fc Sugivamal 

If we allow wi =/= 0, the behaviour becomes even more 
complicated. From Eq. |S1 for wi < 0, lima^oaPx — 0, 
and the dark energy term in Eq. |31 tends to zero. This 
means that all closed cosmologies [flK < 0) will eventu- 
ally recoUapse as da/dt = for some finite a > 1. Open 
cosmologies expand forever, with linia^oada/ dt — \/flK- 
When fix > 0, open models become matter dominated at 



late times, and can go through a phase where fli\i{a) in- 
creases (although obviously pM does not). For wi > and 
fix > 0, the dark energy increases monotonically with the 
scale factor, and it is increasingly unlikely that the matter 
can cause recoUapse before the dark energy takes hold and 
accelerates the expansion of the Universe. The region of 
models which have fix > 0, but that recoUapse therefore 
becomes smaller. 

Solutions that recoUapse are of particular importance 
for the spherical top-hat collapse model: for A cosmolo- 
gies, overdense regions that recoUapse exactly follow the 
behaviour of these models. As we will see later on, for 
general w{a) cosmologies we cannot link spherical pertur- 
bations directly to one of these cosmological models if the 
dark energy does not cluster on the scales of interest. 

Given present observational constraints, the region of 
{flM,flx) space probed in Fig. Qis only of academic in- 
terest, and in the remainder of this paper we focus on the 
subset of models with < flm < 1 and < fix < 1. 

3. linear growth of fluctuations 

Considering the behaviour of homogeneous spherical per- 
turbations provides one of the most simple models for 
the formation of structure in the Universe. The ease 
with which the behaviour can be modelled follows from 
Birkhoff's theorem, which states that a spherically sym- 
metric gravitational field in empty space is st atic and 
is alw ays described by the Schwarzchild metric ijBirkhofj 
Il923l) . This gives that the behaviour of an homogeneous 
sphere of uniform density can itself be modelled using 
the same equations of Section 12 One of the impor- 
tant applications of the spherical perturbation model is 
the derivation of the linear growth rate (pioneered by 
IZel'dovich fc Barenblatt 1958: Peebles 1980. section 10). 
The application proceeds as follows: We consider two 
spheres containing equal amounts of material, one of back- 
ground material with radius a, and one of radius ap with 
a homogeneous change in overdensity. Henceforth quanti- 
ties with a subscript p refer to the perturbation, while no 
subscript relates to the background. The densities within 
the spheres are related to their radii, with 



Ppa 



pa^, 6 = pp/p~l, 



giving, to first order in in 5, 
Up = a(l — S/3). 



(9) 



(10) 



The cosmological equation for both the spherical pertur- 
bation and the background is 



1 d'^a 



Hi 



flMtt-^ + [1 + 3w{a)]flxa^''''^ 



(11) 



where a should be replaced by Up in the matter den- 
sity term for the perturbation. The dark energy density 
px oc a-^*^"', is the same for both the perturbation and the 
background if the dark energy does not cluster. Because 
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Fig. 1. Plots showing the evolution of the matter and vacuum energy densities for a selection of cosmologies with 
different dark energy equation of state (grey lines). For cosmologies with > 0, a = 1 was fixed when the matter 
density and vacuum density were equal. The critical models that border the different types of evolution are shown by 
the black lines. Dashed lines in all plots show the critical models for A cosmologies. The dotted line for wi > shows 
= 0, to emphasis that recoUapse can occur if ^Ix > provided that ^Im >> ^x- 



of this, substituting Eansl^fcllUI into this equation gives, 
to first order in S, 



nna-^ + [3w{a) + ^Qxa^'^"^] } , (13) 



3o u2 -3. 2dad6 
Changing variables from t to a gives 



which can be further simplified to give 



(12) 



da? 



dS 



da 



3 dHnS 
2"*^(") = dh^ 



1 



c? In (5 \ din 5 



d\n( 



dint 



{1 



-[nM{a) + [3w{a) + l]^x{a)] 



(14) 
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Fig. 2. Plots showing contours of constant linear growth 
factor to present day for a selection of cosmologies (grey 
and black lines). The growth factor is normalised to unity 
for = 1, — 0. Increasing w from —4/3 to —2/3, 
decreases the growth factor leading to behaviour more like 
that of an open Universe. Dotted contours are for the fit- 
ting formula of Eq.EOI together with Eans.lT^fcl^ The 
dashed line highlights flat cosmological models. 



This is t he generalisation of Equation B7 in 
IWang fc Ste inhardt (1998) to non-flat cosmologies, 
and is valid for general w{a). 

Eq. 1141 can easily be solved by numerical integration, 
and we show contours of different linear growth factor 



to present day in {flM,^x) space in Fig. [3 for constant 
■w — —2/3, —1, —4/3. The linear growth factor is strongly 
dependent on w, with w > —1 models behaving more 
like open Universes than w < —1 models as the effect 
of the dark energy diminishes. For the linear growth fac- 
tor in a variety of cosmological models with time depen- 
dent w(a) see, for exampl e . iDoran. Schwindt fc: WetterichI 
1I2OOII) : iLinder fc Jenkind (120031) . 

For A co smologies, th e growing mode solution to this 
equation is l|Heathlll977l) 



Dia) = '-^Eia) 



da' 



2'E{a')]3 



(15) 



where E{a) is given by Eq. |3| Although this integral can 
be easily solved numerically, it is comrn o n to use the 
approximation of 'Carroll. Press fc Turner fl992V which 
follows frorn work in .Lightman fc Schechter. (.1990,1 and 
iLahav et all ()l99l|) 



D{a) 



17M(a)^/'-r!A(a) 

^Mia)\ ^ S 1a(q) 
70 



1 



(16) 



A general solution for the growing mode solution in 
dark energy cosmologies, equivalent to Eq. ^1 has yet 
to be found. However, for flat cosmological models, with 
constant w, the solution can be written in terms o f the 
hypergeometric function 2F1 ((Silveira fc Wagalll994^ 



D{a) =a2Fi 



1 w-1 5 

, , 1 , —a 

3w 2w 6w 



M 



M 



Writing the growth index as 

— — = ^Mia), 
a ma 



.(17) 



(18) 



IWang fc Steinhardt ()l998j) use Eq.^for the special case 
of flat cosmologies to give 

3 



5 — w/(l — w) 

3 (l-w)(l-3w/2) 
125 (1 - 6u;/5)3 



l-r!M(a)]. 



(19) 



This ledlBasilakod (l2003l) to ex tend the approximation of 
ICarrofl. Press fc TurneT^1992() given by Eq.[TSlto the 1 

oi w ^ —1 



■ case 



5fiAf(a)a ro / ^a r, f \ 



1 



riM(a) 



1 -1 



il + Anx{a)) 



(20) 



with a given by Eq. [13 and ~ 1.742 3.343w -h 
1.615w^. In Fig. Owe plot the A values required to match 
Eq. 1201 to the true linear growth factor (given by Eq. I17() , 
for flat cosmological models with 0.1 < < 0.9 as a 
function of w (grey lines). The fit of iBasilakosI |2003.) is 
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Fig. 3. Plot showing the true value of A in Eq. 1201 as a 

function of w for 9 flat cosmologies with evenly spread 
between 0.1 and 0.9 (grey lines). For comparison we plot 
the fitting formula of Ba silakoa (,2003,1 (dashed line) and 
for Eq. EH (black line). 



shown by the dashed line. This is a poor fit for w < 
so instead, we propose 

-0.28 



A 



0.08 



0.3, 



(21) 



shown by the black line in Fig. O 

Eq. 1^ has been determined by fitting to flat cosmo- 
logical models with 0.1 < flM < 0.9. For non-flat mod- 
els, the approximation of Eq. |2D1 remains a good fit. In 
Fig. [3 the dotted contours are for the fitting formula of 
Eq. 1201 together with Eqns.El&[2] compared with the 
true value of Dq given by the solid contours. The fitting 
formula fails for Hm << 0.1, but for JIm > 0.1, the maxi- 
mum error (with < rJx < 1) is 3.8% for w = -4/3, 2.6% 
for w = —l and 5.1% for w = —2/3. For c omparison, the 
fitting formula of lCarroll. Press fc Turner! (j^9^ given by 
Eq. ^]is accurate to 2.1% for w — —1 over this range of 

4. the critical density for collapse of spherical 
perturbations 

We now calculate the critical overdensity for collapse of 
homogeneous spherical perturbations at present day in 
a homogeneous dark energy background. The method 
adopt ed is a development of that in IPerc ival ct al.' 
l)2000(l . where the critical overdensity in A cosmolo- 
gies was calculated. S olution scheme s for a n Einstein- 
de Sitter c osmology jG unn fc Gottj Il972|) . for open 
cosmologi es iLacev fc Cola Q993.1 and for flat A cos- 
mologies llEke^'^ ole'^r ^renklll99q) were summarised in 
iKitavama fc Sutol lll996l). See also the solution s cheme for 
A cosmologies used b y iBarrow fc Stein-Schabes (J984) : 
iBarrow fc Saichl iim^ . 

As in Section |21 we consider two spheres containing 
equal amounts of material: one of background material 



with radius a, and one of radius Op with a homogeneous 
change in overdensity. If the dark energy component is 
negligible, as at early times for lima^ow{a) < —1/3, 



/ da 



\d{Hot) 



a 



(22) 



where e is allowed to take any real value. For the back- 
ground, e — Qji = (1 ~ — ^x) is the standard curva- 
ture constant. The matter density Q m is the same for both 
the perturbation and the background as th e two spherical 
region s contain the same mass. Following IPercival et alJ 
l(200nl) . a series solution for a{Hot) in the limit Hgt 
can be obtained given by a = a{Hot)'^^^ + /?(i?oi)'*^^ + 
0[{HQtf/% where 



1/3 



13 = 



3e / 12 



20 \n 



I'M 



1/3 



(23) 



Using the fact that the spheres contain equal mass, 



lim 5{H,t) = -(/?- P,)Wf'^ + 0[{H,tY/\ (24) 
Hot^o a 

Substituting Eq. [521 into Eq. |^ gives that 
lim SiHot) = 



20. 



M 



2/3 



[{\-nM-nx)-ep]{HQtf'^. (25) 



We now extrapolate this limiting behaviour to present 
day, using the linear growth factor to extrapolate the nor- 
malisation of the density field from early times. This gives 



5c = Dq lim 

Hot-*0 



m 



(26) 



where Dq is the linear growth factor to present day. 
Given the normalisation of D{t) adopted in Section 13 
limHot-tQ D{Hot) = a{Hot)'^/^, where a is given by Ea.l^ 
We can therefore write 5r as 



Do 



5n 



M 



[(1 - Om - nx) - ep] 



(27) 



This formula relates the linearly extrapolated overden- 
sity Sc to the initial curvature of the perturbation for cos- 
mological models in which the dark energy becomes neg- 
ligible as Hot — > 0. 

4.1. case 1: cosmological constant 

For A cosmologies, the curvature of the perturbation (and 
energy) is conserved through the Friedmann equation 



da,, 



d{Hot) 



a. 



I I n 2 



(28) 



This equation defines a mini-cosmology, so using the 
methodology of Section[21 the requirement for collapse can 
be seen to be — > 27r2Af^M/4 (compare with Eq. |H1). 
Given a perturbation that collapses, the time taken can 
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Fig. 4. Plots showing contours of constant critical over- 
density for a selection of cosmologies (grey and black lines) 
with constant w. Contours are plotted as a ratio of i^EdS 
as given by Eq. 1^ from 0.96(5EdS to 1.02(5Eds at inter- 
vals of 0.005(5Eds- For the w = —1 cosmology we also plot 
the values determined by accurately solving Eq. using 
standard Elliptical integrals (dotted lines). The difference 
between solid and dotted lines are caused by the numerical 
integration of Eons lSUIl^ll fclS^ The dashed line highlights 
flat cosmological models (as in Fig. [^J. 



be found by integrating Eq. |2H1 which can be reduced to 
a function of standard Elliptical integrals. 



If r^A = 0, Eq. |2H|can be solved analytically and, for 
an Einstein-de Sitter cosmology, the critical overdensity 
for collapse at present day reduces to 



1.686, (29) 
which was first derived bv lQunn (|l972h . 



4.2. case 2: general w{a) cosmologies 

When the dark energy does not cluster, the energy within 
a spherical perturbation is not conserved (a co mprehen sive 
discu ssion of this is given in IWeinberg fc Kamionkowskil 
Although this means that we cannot use Eg. 1281 to 
determine the behaviour of the pertu rbation, we can still 
use th e cosmology equation (Eq. [TT|l . IWang fc Steinhardtl 
l|l998|) provide a boundary value problem for solving this 
second order differential equation, setting the boundary 
at the turn-around time. However, given the relatively 
straightforward initial evolution considered above, it is far 
simpler to set up an initial value problem (in either Op 
or S) t o determine the subsequent behaviou r. In a recent 
paper. Ichiba. Takahashi fc Sueivainal l)2005(l have also de- 
termined the limiting initial conditions for the differential 
equations, working directly from the cosmology equation. 
However, it is perhaps more intuitive to link the critical 
overdensity to the initial curvature of the perturbation 
(Eq.|27I). 

To see how this proceeds from the derivation above, 
suppose that we know Sc- Then Eq. 123 can be used to 
find the initial curvature e^. The initial conditions, Op, S, 
and dap/ dt at some small, but finite a can be obtained by 
substituting into Eans 1221 fc 1251 

Given these initial conditions, the evolution of the per- 
turbation is uniquely specified by three equations: the cos- 
mology equation for the perturbation and the background, 
and the Friedmann equation for the background. For com- 
pleteness we repeat these equations here. 



1 



dh 



ap d{Hotf 
1 d'^a 



a d{Hotf 
da 



^Mttp^ + {3w{a) + l)nxa^''^''A , (30) 



^Ma-^ + i3w{a) + l)nxa^^'''> , (31) 



1 



d{Hot) 



(32) 



By setting the initial value for the solution of these 
equations, we avoid any issues to do with symmetry: for 
general w{a) cosmologies, not only is Eg. 1281 invalid, but 
additionally the evolution of ap in time is no longer sym- 
metric about the point of maximum expansion, because 
the effect of the dark energy is not symmetric about this 
point. This asymmetry is automatically accounted for by 
solving these equations from an initial value. In fact, this 
asymmetry leads to potentially interesting behaviour for 
the spherical perturbations: given the right initial condi- 
tions, perturbations can start to recoUapse (turn-around), 
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but then the repulsive force of the dark energy can cause 
re-expansion l|Chiba. Takahashi fc Sugivaniall2005|) . 

In Fig. ^ we show contours of constant S^, as a func- 
tion of fijv/ and for constant w = — 4/3, —1, —2/3. 
Contours are plotted as a function of Jgds as given by 
Eq. 123 As rix 0, the solutions a symptote towards th e 
open Universe values calculated bv lLacev fc Cold l|l993|) . 
For Qx > 0, we see that decreasing w from —2/3 to 
—4/3 increases the effect of the dark energy. As for the 
linear growth factor, w > —I cosmologies behave more 
like open cosmologies than those with w < — 1: decreasing 
w increases the importance of the dark energy for calcu- 
lating Sc- The critical overdensities for collapse in a va- 
riety of cosmologies with va rying w{a) are compared in 
iMota fc van de Bruckl l)2004|) . 

5. the evolution of the critical overdensity 

The evolution of the critical overdensity for collapse Sc{a), 
is usually defined as follows: for a cosmological model with 
parameters VIm & ^x, <5c(a) gives the overdensity for a 
perturbation that collapses at scale factor a, normalised at 
present day. For example if we had a density field (and as- 
sociated power spectrum) normalised at present day, then 
dc(a) (where a does not necessarily equal 1) relates to 
spherical perturbations in this density field that collapse 
at scale factor a. Sc is the particular case for perturbations 
that collapse at present day: perturbations that collapse 
earlier obviously have to be significantly more overdense. 
To calculate Sc{a), when integrating the Friedmann equa- 
tion fEa.l28|l for A cosmologies, or numerically integrating 
Eans. lSnil31l fcl^we would look for collapse at scale factor 
a rather than a = 1. 

The rather weak evolution of the critical overdensity 
as a function of cosmological model (Fig. 0)) , means that 
the linear growth factor can be used to approximate Sc{a)'. 
If Sc is constant along a particular cosmological track then 
the evolution of Sc{a) is purely driven by D{a)^^ . The only 
change in 5c{a) between two collapse times is caused by 
the change in overall normalisation of the field, which can 
be seen by considering Eq. ESI for two cosmologies along 
a single cosmological track. The most obvious choice for 
the normalisation is (5EdSi so the approximation will be 
correct in the limit as a ^ 0, 



<5c(a) 



_Do_ 

D{a) 



JEdS- 



(33) 



In order to demonstrate this approximation, in Fig. [3 
we plot the ratio between the true critical overdensity 
Sc{a) to the approximation given by Ea. 1331 for three cos- 
mological models with fijvf = 0.3, ^x = 0.7, but with 
constant w = -2/3,-1,-4/3. As a — > 0, the approxima- 
tion becomes more accurate. At a = 1 (redshift 0), the 
ratio is the same as the value of the surface contoured in 
Fig.0] We also plot the approximate value of 5c{a) calcu- 
lated from Eq.lJH but with D{a) and Z?o calculated using 
the approximation given by Eq.lJOl Eq.|^& Eg. 119! fgrev 
lines) . The error from using the approximation of Eq. 133! 




Fig. 5. Plot showing the ratio between the critical over- 
density for collapse Sc{a) and the approximation given by 
Eq. 123 as a function of redshift for three different cosmo- 
logical models (black hues). VLm = 0.3 and Vlx = 0.7 
are assumed, while w is assumed to be constant with 
w — -2/3,-1, or —4/3. The grey lines shows the ra- 
tio between the true critical overdensity for collapse and 
the approximation, where the linear growth factor has it- 
self been approximated by Eq. [201 The error in using the 
Einstein-de Sitter critical overdensity is of the same or- 
der as the error in the approximation of using the fitting 
formula of Eq. ^ 



is of the same order as that from assuming the approxi- 
mation to the linear growth factor for these cosmologies. 



6. the mass function 

The usefulness of the spherical model was emphasised 
when Press & Schechter considered smoothing the ini- 
tial density field to determine the r elative abundances 
of pe rturbations on different scales ijPress fc Schechter! 
Il974t PS). Whe n combined with the critical overdensity 
for collapse this provided a statistical model for the for- 
mation of structure in the Universe: smoothing the fluc- 
tuations leads to the masses of collapsed objects, while 
the spherical perturbation model gives the epoch of col- 
lapse for those perturbations that are sufficiently dense. 
Obviously such a simple model will fail in detail, par- 
ticularly given the known complexities of asymmetrical 
gravitational collapse, and num erical simulations have 
now quantified these problems (|Sheth fc Tormenlh999l: 
{jenkins_ct al. 2001). However PS theory has been incred- 
ibly successful and arguably still provides key insight into 
the processes at work in structure formation. 

Two alternative formalisms are often considered for 
the collapse of perturbations in PS theory 

1. The overdensity field is assumed to grow with the lin- 
ear growth factor, and when perturbations reach the 
critical overdensity they are said to have collapsed. 
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12 14 
logj(,(mass / Mq) 

Fig. 6. Plot showing the predict ed mass function calcu- 
lated using the fitting formula of lSheth fc Tormenl l)l999(l 
calculated for Qm = 0.3, fix = 0.7 for three differ- 
ent values of w, and at three epochs corresponding to 
a = 1/3, 1/2, 1. Because the power spectrum is normalised 
at present day, and 6^ is only weakly dependent on cosmol- 
ogy, then there is little difference between the predicted 
mass functions for a — 1. As we go further back in time 
the difference becomes more severe because of the differing 
linear growth factors. 

2. Each overdense region is considered to be spherical and 
its collapse time is calculated as in Section 

Given the discussion in Section|Sl it is easy to see that the 
first formalism, corresponding to a growing field, matches 
the approximation to Sc{a) using the linear growth factor 
given by Eq. 1331 The second formalism, which is adopted 
in the following analysis, corresponds to using the correct 
6r. ( a) for the spheri cal model. 

iMa et all l|l999l) considered the effect of quintessence 
on the mass transfer function. They provided fitting for- 
mulae for the ratio between the quintessence and A 
cosmologies. However, if the dark energy only clusters 
on very large scales, the transfer function is only al- 
tered on these scales. If the power spectrum is nor- 
malised to tTg (the rms density fluctuation on scales of 
8 h Mp c~^), then the scales usually of interest are not af- 
fected l)Lokas. Bode &: HoffmanlEool . 

To demonstrate the effect of the dark energy equa- 
tion of state on the mass function, we plot the cumu- 
lative mass func tion A^(> M) . calcu lated using the nu- 
merical fit of .Sheth fc TormCT l)l999(l for 3 different cos- 
mologies and 3 different epochs in Fig. |H1 To determine 
the power spectrum we have used the fitting formulae of 
lEisenstein fc Hul l)l998|) with = 0.3, nt/flM = 0.15, 
(Tg = 0.9, ris = 1. The critical overdensity for collapse 
for w = —2/3,-1 & —4/3 and fix = 0.7 was then cal- 
culated for a = 1/3,1/2,1, corresponding to redshifts 
z — 2,1,0. As expected, because the critical overdensity 
for collapse is only weakly dependent on cosmological pa- 
rameters (Fig.0}, at a = 1 (the epoch at which the power 



spectrum is normalised) we see very little difference in 
the predicted mass functions for different cosmologies. If 
the normalisation of the power spectrum had been con- 
strained at a different epoch (for example by CMB fluctu- 
ations), then this would not be correct. The evolution of 
the mass function is strongly dependent on w because of 
the effect on the evolution of (5c (a) (Section^ through the 
linear growth factor (Fig. ^ . Consequently, determining 
the mass function at redshifts other than that used to nor- 
malise the power spectrum offers a stronger possibility of 
measuring w{a). We will discuss the evolution of structure 
growth further in the next Section. 

In order to compare with the observed cluster counts, 
we must convert from comoving position to observed an- 
gular position and redshift. The number of sources with 
mass > M per unit solid angle in a redshift slice dz is 
given by 



N{> M,Az) 



z+Az 



dzN{>M)d%^^ 



dd-n 



dz 



where the proper distance dprop is given by 
c da 



(34) 



(35) 



This correction to the mass function reduces the signifi- 
cance of w at low redshifts for masses of order 10^** 
()Solevi et al...2005,) . 

7. rate of structure growth 

As discussed in the previous Section, the present day mass 
function does not provide a good test of w, although its 
evolution in time does. The Pre ss-Schechter mass func- 
tion (and the fitting formula of ISheth fc Tormenl 1 19991) 
is independent of epoch when written as a function of 
V — 5c/ iM, where um is the rms fiuctuation in the ini- 
tial overdensity field on a scale corresponding to mass 
M . Consequently, it is easy to see that the evolution of 
the mass function is dependent on the rate at which 5c 
changes, d5c{a) / da. It is this quantity that we are testing 
by comparing mass functions at different epochs. 

Extended Press-Schechter theory uses the smoothed 
initial overdensity field t o create an analy t ic model for 
the build- up of structure ('Bon d et al]ll99 it lBowe r"l99l': 
iLacev fc'~Cole,199 3'). in addition to providing a statistical 
model for the mass function. Within this extension to the 
theory, the rate at which structures grow i s again driven 
bv d5Ja)lda JPerciyal fc Milleil Il999l: IPercival et alJ 
ES^ ICo1e et a,1.l l2nnnt iMiHer. Percival fc CroomI \2m^. 
Appendix A of iMiller. Perc ival fc CroomI l(20fl5l) showed 
that the approximation of Eq. 1331 can be easily extended 
to d5c{a)/da for A cosmologies. 

Here, we extend this analysis to consider d5c{a)/da 
for w ^ —1 cosmologies. In Fig. |7| we show d5c{a)/da 
for three flat cosmological models with VIm ~ 0.3, and 
w — —4/3, —1, or —2/3 calculated using the derivation of 
Section ^21 (black fines) and the fit of Eq.jSSKgrey lines). 
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Fig. 7. Plot showing the derivative of the critical overden- 
sity for collapse with respect to the scale factor d5c{a)/da 
(black lines) and the approximation calculated from Eg. 1331 
(grey lines), as a function of redshift for three different cos- 
mological models. ^Im = 0.3 and VLx = 0.7 are assumed, 
while w = —4/3,-1, or —2/3. The approximation is so 
good that the grey lines are almost completely hidden be- 
hind the black lines. 

In fact, the fit is so good that the grey lines are hardly 
visible in this plot. As expected d5c{a)/da is strongly de- 
pendent on the dark energy equation of state. 

8. virialisation 

The inhomogeneous nature of true perturbations means 
that they do not collapse to singularities, but instead sta- 
bilise at finite size. For A cosmologies it is possible to use 
energy considerations to determine the final radius and 
density of the virialised perturbation 



ijLahav et alJll99Ll . 

The extension of this analysis to more general dark en- 
ergy c osmologies h as been considered by a number of au- 
thors llWang fc Stein hardt 1998; HorcUou & Bcrec 200^ 
lMa,or fc T.a,ha,vll2nn5t). 

Within a perturbation, the potential energy due 
to the matter Ug and dark energy Ux ar e given by 
llHorellou fc Bergell200,'il: iMaor fc Lahavll200,'j) 



Ug 



5R ' 



AwGM 



10 



(36) 



These can be calculated from the Poisson Equation with 
pressure term. Note that there is some confusion in the 
literature about the exact form of Ux and the [1 -I- 3w{a)] 
term has sometimes been neglecte d in the p ast, although 
it is included in more re cent work ijBattve fc W cUer 2003: 
iHorellou fc Bergell200,'^^ . Even without this term, the dis- 
cussion below about the lack of energy conservation within 
the perturbation remains valid, although the numerical re- 
sults will obviously change. 

For dark energy cosmologies the virial theorem, that 
a system with potential energy t7 cx i?P virialises with 



temperature T = pU/2, holds and T — —^Uq + Ux at 
the epoch of virialisation. Assuming conservation of total 
energy between turn-around at ata and virialisation Cvir 
gives 

1 



C/G(ata) + Ux{at^) - 2^G(avir) + 2Ux{a.ir). (37) 

Sub stituting the defini t ions given by Eg. 15^ leads 
to llWang fc Steinhardtl Il998t iHorellou fc BergJ 120051: 
iMaor fc Lahav 2005') 



( 



l-[l + 3u;(ata)]|ja;+[l 



3w(av 



qx^ _ 1 

y 2 



(38) 



where 



Rta 



[1 -I- (5(ata)]f^Af (ata) ' 



y 



pxjata) _ a. 
Px(avir) 



/(ata) 



/(Qvi, 



(39) 



Here q gives the ratio of the dark energy density to the 
matter density in the perturbation at turn-around. For 



w 



-1, 



1 



. and Eq. [ 

iLahav et all l)l99j) . For w 
X = 1/2. 

The above derivation 
tion that total energy is 
tion between turn-around 



reduces to the formula of 
— 1/3 or ilx = 0, we find 



was based on the assump- 
conserved for the perturba- 
and virialisation. This is a 
common assumption in the literature for cosmologies 
in which the dark energy rem ains homogeneous on the 
scales of the perturbations (e.g. Wang fc Steinhardt"l998l 
Weinberg fc Kam ionkowski 2003; Battyc fc WcUcr 2003; 
iHorellou fc Bergel l2005). Even in work that considers pos- 
sible dark energy clustering, energy conservation has pre- 
viously been assumed in the hom ogeneous case (see th e 
discussion leading to equation 26 of lMaor fc Lahavll2005|) . 
However, the potential energy of the matter due to the 
presence of the dark energy Ux {U12 in the notation of 
iMaor fc Lahavl2005|) is dependent on the dark energy den- 
sity in the perturbation px (a) . The evolution of this den- 
sity lies with the background, rather than the perturba- 
tion, so the total energy is not expected to be conserved: 
this is why we could not write down a Friedmann equa- 
tion for the perturbation in Section ^ Obviously, if the 
total energy of the perturbation between turn-around and 
virialisation is not conserved, then Equation 1371 does not 
hold. 

In fact, we can consider two extreme situations. In the 
first, as above, the total energy at Ovir is the same as at 
flta- In the second, we assume that the perturbation did 
not change size between the two epochs. Following the 
second assumption, the total energy in the system at Ovir 
would have been altered from that at Ota by the change 
in the dark energy potential Ux,i — Ux.2, where 



Ux,i = [1 -t- 3u;(ata)] 



AttGM 
10 



(ata)-Rta) 



(40) 
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Fig. 8. The total energy in an ideal homogeneous spher- 
ical perturbation (as modelled in Section 21 , undergoing 
collapse to a singularity at a = 1, relative to the energy 
at turn-around (black lines). CIm — 0.3 and fix = 0.7 
are assumed, with constant w at the values shown. For 
general w, the total energy is not conserved. If the pertur- 
bations had remained at the turn-around size throughout 
their evolution, then the total energies would have evolved 
along the grey lines. The actual evolution in total energy 
lies between the line of constant total energy and the line 
of constant size for each cosmology. 



Ux.: 



[1 + 3w{a^ 



, . 4ttGM , 



)R 



(41) 



We call using the dark energy potential Ux,i to calculate 
the total energy in the perturbation at virialisation "fixing 
the energy at turn-around" , and using Ux,2 "fixing the 
energy at virialisation" . This is demonstrated in Fig. [SJ 
where we plot the total energy in an ideal homogeneous 
perturbation undergoing collapse to a singularity at a = 1 
relative to the energy at turn- around (black lines), for a 
number of cosmologies. The total energy was calculated by 
combining the analysis of Section^with the definitions of 
potential energy of Eq. |2H Obviously if the total energy 
were constant we would find a horizontal line, as is the 
case when w = —1. If the perturbation did not evolve 
in size, with the radius fixed at the turn-around value, 
then the evolution of the total energy is shown by the 
grey lines. The actual total energy within the perturbation 
lies between these two extreme cases for each value of w 
assumed. 

For a more realistic perturbation that underwent viri- 
alisation rather than collapsing to a singularity, provided 
that the dark energy potential varies monotonically be- 
tween Hta and flvir, then the change in total energy in the 
system between turn-around and virialisation should still 
lie between these two cases: some of the change in total en- 
ergy due to the dark energy remaining homogeneous will 
be converted into kinetic energy, and some will be "lost" 
to the background. 



energy fixed at turn-around 
energy fixed at virialisation 




-1.5 



-0.5 



Fig. 9. Upper panel: the ratio between the virialisation 
radius and maximum turn-around radius for a spherical 
perturbation that virialises at a = 1 (black lines) or a = 
1/2 (grey lines), assumed to be the same as the epoch 
predicted for collapse of a homogeneous perturbation. We 
assume that VLm = 0.3, Vtx = 0.7, and plot results as a 
function of constant w. The dashed and solid lines show 
two different calculations assuming conservation of energy 
based on the dark energy density either at turn-around, 
or at virialisation. Lower panel: as upper panel, but now 
showing the overdensity of the virialised system. 

The second scenario, using Ux,2 to "fix" the energy at 
virialisation alters Eq. |2H1 to give 

1 - [1 + 3w;(avh.)]|-) a; + [1 + 3u.(avir)]— = ^, (42) 
22// y 2 

with X, y, Sz q defined as before. For w = —1 t his eq uation 
again reduces to the formula of iLahav et alJ l|l99lh . 

In the upper panel of Fig. Owe plot i?vir/-Rta in these 
two cases in a cosmology with ^Im ~ 0.3, fix — 0.7, as a 
function of constant w. We assume that virialisation oc- 
curs at the collapse epoch, set to be either the present 
day or a — 1/2, for the homogeneous spherical model. For 
constant u; = — 1/3 we find Rvii/Rta = 1/2 as expected. 
For cons tant w = —1 , the t wo solutions converge to the 
result of iLahav et"all l|l99l|) . For -1 < w < -1/3, the 
component of the energy of the perturbation in the form 
of dark energy changes from turn-around to virialisation, 
and both energy conservation arguments predict different 
values of i?vir/-Rta, although the trend with w is the same 
for both. However, for w < —1, the solutions diverge. It is 
easy to see why: for w << —1, the dark energy becomes 
increasingly important at late times in the evolution of the 
Universe. For a solution that collapses at present day (for 
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example), turn-around must happen at an epoch where 
the dark energy has yet to be cosmologically important. 
It is only between turn-around and virialisation that the 
dark energy becomes important. Setting the energy to be 
constant based on the dark energy density at virialisa- 
tion simply continues the trend from —1/3 > w > — 1, 
with px increasing as w decreases, leading to a decrease 
in the potential due to the dark energy Ux, and a smaller 
i?vir/-Rta ratio. However, setting the dark energy density 
at turn-around has the opposite effect: px decreases, lead- 
ing to Ux and i?vir/-Rta increasing. For perturbations that 
turn-around, but do not collapse, the dark energy density 
increases so rapidly between turn-around and virialisation 
that the predicted i?vir/^ta keeps increasing (although 
virialisation is never reached because the perturbation can 
never stabilise). 

In the lower panel of Fig. O we plot the overdensity at 
virialisation Avir, calculated for the two values of the total 
energy. As can be seen, the result has a very similar depen- 
dence on w as described above for i?vir/-Rta. It is clear that 
further analysis is required before measurements that de- 
pend on Avir can be used to constrain w. The answer may 
lie in numerical simulations of the time evolution of the 
pe rturbation leading to virialisation, po s sibly of the form 
of lEngineer. Kanekar fc Padmanalph ^n] j2QQ(j |). or possi- 
bly standard N-body sim ulations l|Meneghetti et alJl200,'Tt 
iBartelmann et al.ll2005|) . Such simulations will need to 
consider the evolution of total energy with time as well 
as the changing epoch of virialisation: given the possible 
strong effect of the dark energy, it seems clear that the as- 
sumption that virialisation of the perturbation occurs at 
the collapse time predicted in the homogeneous case will 
also break down in addition to the conservation of total 
energy. 

9. conclusions 

The presence of dark energy alters the way in which cos- 
mological structures grow, thus providing an observational 
signature that is complementary to geometrical effects. 
The structure growth is dependent on the normalisation 
of the dark energy density Qx, its equation of state, and 
sound speed, which determines how this component clus- 
ters. In this paper we have only considered the effect of 
the first two of these properties, assuming that the dark 
energy remains homogeneous on the scales of interest due 
to a high sound speed. The spherical top-hat model has 
been used to determine the linear growth rate and non- 
linear collapse overdensity threshold. The equations pro- 
vided have allowed for general ■w{a). We have also con- 
sidered the statistics of observed structures through the 
mass function and its evolution. Finally, the virialisation 
of perturbations has been considered and a new argument 
has been presented demonstrating the importance of the 
lack of energy conservation within a perturbation. It is 
clear that more work is required before observations that 
use arguments based on the energy in perturbations can be 
used to constrain the dark energy. However, there is clearly 



tremendous potential for future observations to detect the 
cosmological effects of dark energy in sufficient detail to 
pin down its properties. Proving that w ^ —1, or indeed 
that w = —1 would be an exciting result, and remains one 
of the goals of modern cosmology. 
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